(N 



Cooperative localization using angle of arrival 
measurements: sequential algorithms and 
non-line-of-sight suppression 

Bharath Ananthasubramaniam and Upamanyu Madhow 
Electrical and Computer Engineering Department 
University of California, Santa Barbara 
Santa Barbara, CA 93106 
E-mail: bharath@engr.ucsb.edu, madhow@ece.ucsb.edu 

Abstract 

. We investigate localization of a source based on angle of arrival (AoA) measurements made at a geographically 

dispersed network of cooperating receivers. The goal is to efficiently compute accurate estimates despite outliers in 
the AoA measurements due to multipath reflections in non-line-of-sight (NLOS) environments. Maximal likelihood 
(ML) location estimation in such a setting requires exhaustive testing of estimates from all possible subsets of 
"good" measurements, which has exponential complexity in the number of measurements. We provide a randomized 
algorithm that approaches ML performance with linear complexity in the number of measurements. The building 
block for this algorithm is a low-complexity sequential algorithm for updating the source location estimates under 
line-of-sight (LOS) environments. Our Bayesian framework can exploit the ability to resolve multiple paths in 
wideband systems to provide significant performance gains over narrowband systems in NLOS environments, and 
c/2 ' easily extends to accommodate additional information such as range measurements and prior information about 

location. 
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\q \ I. Introduction 

^ I In this paper, we investigate cooperative localization of a source transmitting a known signal using a network 
of geographically dispersed receivers (detectors) using angle of arrival (AoA) measurements. We consider AoA 
measurements, since they only require that each receiver has a calibrated antenna array with known orientation and 
location, and that the receivers can coordinate to pool all the AoA estimates corresponding to a given source (e.g., 
using coarse timing synchronization between the receivers to associate AoAs for the signal from a given source at a 
cd ■ given point of time). While the geometry of AoA-based source localization is straightforward for line-of-sight (LOS) 
environments, our primary goal in this paper is to develop efficient localization algorithms for non-line-of-sight 
(NLOS) environments. 

While the problem considered here is of broad applicability, our primary motivation is localization in sensor 
networks, where a network of collector nodes (i.e., data gathering nodes with some advanced capabilities such as 
AoA estimation) collaboratively decode and localize transmissions from elementary microsensors that have data 
to send. Since sensors communicate only when they observe an "interesting" event, without prior coordination 
with the collector nodes, this sensor-driven paradigm (proposed in [1 ]), allows drastic reduction in the microsensor 
functionality and communication energy costs. By locating the sensors transmitting in response to the event, the 
location of the event can be estimated. Our goal, therefore, is to locate the transmitting sensor in NLOS propagation 
conditions. Note that the term 'sensor' is not a reference to receivers, as is sometimes the usage in the source 
localization and array processing literature. To avoid confusion, we use the term "source" to refer to the source of 
the transmitted signal, and the term "receiver" for a device that receives this signal and produces an AoA estimate. 
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This problem is also applicable to asset tracking using active radio frequency identification (RFID) tags, in which 
tags periodically or intermittently transmit a signal to be used to localize them |2), 0, |4). 

In the preceding scenario, other measurements such as time-difference-of-arrival (TDOA) [2] and received signal 
strength (RSS) [5] could also be used for localization. TDOA measurements require stringent synchronization 
between the receivers, while RSS measurements are highly sensitive to the propagation environment, and typically 
require extensive calibration. On the contrary, AoA measurements offer the prospect of providing accurate localiza- 
tion without tight coordination among receivers or extensive calibration, at least in LOS environments. In this paper, 
we investigate whether this promise can be realized in more realistic NLOS environments, with the understanding 
that AoA-based performance can be augmented by incorporating information from TDOA and RSS measurements, 
if available. 

We previously proposed one possible receiver design (including timing acquisition and array processing algo- 
rithms) to generate these AoA estimates (T). In this paper, we abstract away such details in order to focus on 
the fundamentals of AoA-based localization. To this end, we propose statistical models for AoA measurements 
that model the performance of a broad class of AoA estimators in the literature under a variety of propagation 
environments. These models permit performance comparisons between our proposed algorithms and fundamental 
limits such as the Cramer-Rao Lower Bound (CRLB). Our main contributions are as follows: 

1) As a building block for localization in NLOS environments, we develop a scalable sequential algorithm 
for aggregation of AoA estimates from multiple receivers to estimate the source location in LOS scenarios 
(minimal multipath scattering). The receivers only need to exchange source location and covariance estimates, 
where the source location estimate is updated with the local AoA measurement before it is forwarded to the 
next receiver. This approximately maximum likelihood (ML) algorithm has linear complexity in the number 
of receivers. We use the CRLB to estimate the location uncertainty as a function of the coverage area, the 
variance of the AoA measurements, and the number of receivers. 

2) We model the AoA estimates due to multipath components with AoA measurements far away from the 
LOS path as outliers. Assuming that at least some of the receivers see near-LOS paths, ML localization 
requires brute force elimination of outliers by considering all possible subsets of AoA measurements, which is 
excessively complex. We propose a randomized algorithm with 0(MN 2 ) complexity for outlier suppression, 
which employs M randomly initialized instances of the sequential algorithm. Here M is chosen to achieve 
a given probability of localization failure, and does not grow with N, the number of receivers. 

3) The proposed algorithms are numerically shown to achieve the CRLB and the ML performance. We quantify 
the performance penalty due to the presence of outliers, and provide examples that illustrate the performance 
advantage of wideband systems, which can resolve LOS and NLOS paths, relative to narrowband systems 
without the capability for providing such resolution. Also, our Bayesian framework permits integration of 
other sensing modalities such as RSS-based range measurements for localization in a NLOS environment. 

An outline of these preceding results was presented previously at a workshop [6]. This paper significantly expands 
upon that work, including detailed derivations of the algorithms and performance bounds, as well as a more 
comprehensive set of simulation results. 

Related Work: There is a rich history of literature on source localization using a variety of measurements, 
including TDOA, AoA and RSS that depend on LOS channels between the source and the receivers with no 
multipath and consider errors due to measurement noise alone. A detailed discussion of the body of standard 
localization techniques using multiple modalities, their limitations and practical considerations are presented in Q, 

n. 

Our approach to localization in NLOS environments draws upon the significant literature on handling outliers 
J9). One approach is to identify and remove outliers before location estimation. In ifTOll . a generalized likelihood 
ratio test is used to identify NLOS estimates; this assumes at least partial statistical knowledge about the NLOS 
components and does not utilize the inherent consistency between measurements corresponding to a single source. 
Prior knowledge of the NLOS characteristics is used to eliminate outliers in TDOA and AoA measurements in 
iTTTI . but requires exploring all subsets of measurements, which is not scalable. The skewed distribution of range 
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or time of flight (TOF) estimates due to a positive bias introduced by NLOS propagation is used in |[T2l to identify 
and correct for NLOS errors. Correction of NLOS time of arrival estimates prior to localization is accomplished by 
Kalman filtering in 1 1 3 ] and soft combining of pseudo-range measurements in [14]. The residuals from locating the 
source using a least-squares algorithm on a set of AoA measurements are compared against a threshold to identify 
and eliminate outliers in |[T5l . More recently, Yu et al. |[T6l formulate a Neyman-Pearson test to only identify 
the outliers in different combinations of localization measurements such as AoA, ToA and RSS-based ranging 
and similarly, Giivenc et al. [17 | identify NLOS components in multipath ultrawideband (UWB) channel metrics 
using joint likelihood ratio tests. Al-Jazzar et al. |[18 | design a nonlinear optimization algorithm with nonlinear 
constraints to compute the locations of the unknown scatterers to identify the NLOS components. Giivenc et al. 
present a detailed survey of UWB TOA-based NLOS mitigation techniques in |fT9ll . 

In this paper, we adopt the alternative approach of robust location estimation while simultaneously limiting the 
effect of outliers (called NLOS "mitigation"), thus eliminating outliers "on the go." To the best of our knowledge this 
approach has only once been used with AoA measurements in ll20l and similar work for range/TOF measurements 
includes (2T|, 11221 . ll23l . Tang et al. ll20l formulate the localization using AoA and ToA as an optimization with a 
geometrically-constrained objective function that has a higher computational complexity of at least 0(N 3 ). More 
significantly, they only consider angular spreads due to a disk of scatterers around the transmitter that results in a 
simple Gaussian model of AoA spread. Thus, this AoA distribution is equivalent to having Gaussian errors in the 
LOS estimate and does not account for blockage of LOS. Chen ||2~T1 proposes a weighted least-squares localization 
using range/TOF measurements, where the weights are iteratively varied to assign the lowest weights to outliers and 
highest to the LOS estimates. The main drawback of this approach is the 0(2^) complexity due to a combinatorial 
exploration of subsets of "good" measurements, which does not scale. Venkatesh and Buehrer 11221 present a linear 
programming approach (of approximately 0(N 3 ) complexity) to mitigate NLOS using ToA estimates in a UWB 
system that utilizes the positive bias in NLOS range measurements. 

Our approach shares many similarities with Casas et al. ||23l . Casas et al. perform closed-form trilateration of 
random triplets of measurements (similar to the M randomizations in our method) followed by a selection of 
the location estimate with the smallest median residual, which has an O(MN) complexity. They propose a slight 
improvement by performing a multilateration after identification of the outliers from the previous step, incurring an 
0(N 2 ) complexity. However, localization even without outliers using AoA measurements does not have a closed 
form solution in terms of more than two noisy AoA measurements in two-dimensions (and more than three noisy 
AoA measurements in three-dimensions), i.e., in an over specified system, due to the trigonometric dependence 
of AoA on the source location. This nonlinear dependence of AoA on the source location makes the algorithmic 
complexity inherently larger than for the time of flight measurements considered in ll23l . Our sequential algorithm 
offers a low-complexity approach to localization, which can be viewed as a generalization of ||23l : we naturally 
grow "good" subsets of measurements to include many LOS measurements, without needing a separate phase 
for localization using inliers. Furthermore, our Bayesian framework permits incorporation of other localization 
modalities, such as RSS-based range measurements, with little change to the algorithm (see Figure H). 

The rest of the paper is organized as follows. AoA measurement models under different propagation environments 
are described in Section |II] The sequential algorithm for localization in LOS scenarios, along with performance 
benchmarks, are presented in Section |III] We derive an extension to the sequential localization to suppress NLOS 
AoA estimates in Section ITVl which is motivated by the structure of the ML estimate. The proposed algorithms 
are numerically investigated in Section fVl Section |VT] contains concluding remarks. 

II. Models of AoA Estimates 

Each receiver estimates the direction of arrival of the source transmission using an antenna array. In this section, 
we discuss models for the estimated AoA that capture the effects of LOS and NLOS propagation environments. 
These models are then used to design source localization algorithms in Sections HIT] and llV] For simplicity, we 
assume a linear array, although our framework generalizes to arbitrary antenna arrays, as long as the array manifold 
is known. 
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Classical AoA estimation techniques such as MUSIC and ESPRIT [24] are designed to separate uncorrected point 
sources assuming LOS propagation from the source to the receivers. This point-source propagation model does not 
account for multipath scattering encountered in many deployment environments. Multipath adds highly correlated 
AoA components to the received signal, and extensions to deal with such effects using techniques such as spatial 
smoothing have been proposed at the cost of lower AoA resolutions and poorer source separation capabilities. 

Here, we adopt the approach taken in [25 1 , 1126*1 , ||27l , lf28l , lf29l , where a series of models for AoA estimation 
in the presence of scattering have been proposed. The propagation is characterized by a mean arrival angle 
(corresponding to the true bearing) and a spatial spreading parameter (quantifying the spatial uncertainty caused 
by multipath). The resulting AoA at the antenna array is a random variable, typically modeled as Gaussian lf28l or 
Laplacian [30]. Drawing on these ideas, we propose the following models: 

LOS model: We characterize the spatial spreading in LOS propagation scenarios by zero-mean symmetric finite 
variance "noise" models, such as the Gaussian and Laplacian. AoA estimation under LOS in the presence of additive 
white Gaussian noise results in zero-mean Gaussian errors, whose variance depends on the signal-to-noise ratio OT1 . 
Additional errors in the AoA measurements due to local scattering in the vicinity of the source are also modeled 
using zero-mean symmetric distributions. Therefore, Gaussian and Laplacian error models seamlessly transition 
between scenarios with and without local scattering for different values of the spatial spreading parameter. These 
models represent situations where the received signal has a strong LOS component, together with a limited amount 
of scattering. For a source at location X_ along a true bearing 0{X), the Gaussian LOS model with local scattering 
is 

7T 7T~ 
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where a 2 represents the spatial extent of scattering, Q(t) = J°° exp(— t 2 /2)dt/ ^/2tt is the normal tail distribution 
and all angles are measured with respect to the antenna broadside. The permissible angles are within ±7r/2 of 
the antenna broadside and the Gaussian density is truncated to reflect this. As the spatial spread a 2 increases, 
the density become progressively long tailed and tends toward a uniform density. The Laplacian LOS model with 
spatial spreading factor a has heavier tails than the Gaussian model: 



PLap{ e/x) = — I ex P L^em) , e e 

fLap \ / j 2 ^ ex P (-|7r/(V2a)|) \ o J 
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LOS blockage model: Environments where the LOS path to a receiver is blocked by structures, such as building, 
hills or trees, are also relevant. As a result, the received signal is composed exclusively of multipath components 
that deviate "far" from the LOS path. Therefore, AoA measurements made from these scattered and reflected paths 
alone are fairly uncorrelated with the true bearing of the source. An AoA estimate drawn uniformly from the 
feasible set aptly models such a scenario: 



Pnlos(#/2Q = -, 9 £ 
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When there are significant contributions from both the LOS path and the multipath components from other 
directions, the model for the AoA estimates depends on the capability of the receiver to resolve these contributions 
spatially. If the source signal has a large enough bandwidth, the multipath components can be resolved in time, and 
separate AoA estimates can be obtained for each path. If the bandwidth is insufficient to resolve the paths, then 
the receiver may only be able to obtain a single AoA estimate. We model these two scenarios separately. 

Narrowband multipath model: With narrowband source transmissions, the receiver is capable of resolving the 
arriving combination of LOS path and NLOS multipath in the spatial domain only. For receivers with relatively 
small number of antenna elements, the receiver can only measure the AoA of the strongest arriving path (or 
strongest superposition of paths). Hence, each receiver produces a single AoA estimate that, depending on the 
relative strengths of the LOS and NLOS components, will provide an estimate close to or very "far" from the true 
bearing of the source (outliers). Accordingly, we model a typical narrowband scenario as follows: Let a be the 
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fraction of receivers that experience reflected and scattered multipath significantly stronger than the LOS path. In 
this event, the AoA appears to correspond to the LOS path being blocked and hence, is drawn from the worst-case 
NLOS model in Q. In the remaining instances, when the LOS path is strong, the AoA is drawn from one of the 
LOS models, for instance CO. Thus, we arrive at the following narrowband model: 



Pnarrowband(#/^) = a PNLOS (0 / 2L) + (1 - a) p G aus(# / 'K) ■ 



(4) 



This model also represents the least favorable multipath environment, where the LOS component is either present 
and significant, or is blocked (completely absent), and serves as our candidate model for developing the NLOS 
suppression algorithm. Note that distributions with heavy-tails, such as the Laplacian (in (O) or a Cauchy model, 
can also be used to represent the narrowband setting with outliers. We investigate these alternate models in Section 

EE 

Wideband multipath model: With wideband source transmissions, the receiver has an additionally resolve paths 
in time as well. When the scattered or reflected paths are sufficiently temporally or spatially separated from the 
LOS path, the receiver is capable of resolving the LOS path, and possibly multiple NLOS paths. We model the 
multiple AoA estimates produced by the receiver as follows: the AoA estimate corresponding to the LOS path 
is represented by an LOS model such as CO or (|2), while the remaining estimates are drawn uniformly from the 
feasible set, as in the worst-case model with LOS blockage. Thus, if a receiver resolves L paths, the resulting AoA 
estimates are generated according to the following distribution: 



Pwideband^ 1 ), 6™ . . . / X) = PGsas (0^/X) PsuasV® / 20 • • • PNLOS (9&/X) . 



(5) 



Assuming that the LOS path is not blocked, it is resolved by a receiver in the wideband scenario. However, it is 
not known a priori which of the L AoA estimates corresponds to the LOS path. If the LOS path is blocked, then 
none of the L resolved AoA estimates might be close to the true source bearing. 
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cartesian coordinates 
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Aggregated and remaining receivers 
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Fig. 1. (Left panel) One possible sensor network layout with a circular field and 8 receivers on its perimeter: The geometry used to compute 
the "bootstrap" source estimate from AoA estimates of receivers Rl and R2 is shown. (Center Panel) The flowchart of the sequential 
localization algorithm using AoA measurements (with no outliers) presented in Section IIII-BI (Right Panel) The flowchart of the outlier 
suppresion algorithm in Section H V-B I that uses the sequential algorithm in the Center Panel. 



III. Localization in LOS scenarios 

In this section, we consider scenarios with only LOS propagation, i.e., no LOS blockage, where the spread in 
the AoA estimates is caused only by local scattering. We present an algorithm for sequential aggregation of the 
available AoA estimates (generated according to the LOS model in Section [TT]) to produce the estimate of the 
source location. Each receiver performs linear minimum mean squared error (LMMSE) updates on the prior source 
location estimate (SLE) (received from a previous receiver) using its own AoA estimate, before passing the updated 
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estimate to the next receiver. This process is continued until all the available AoA estimates is aggregated into the 
SLE. 



A. LMMSE Updates and Coordinate Transformations 

At each receiver, we desire linear updates (to keep the computational burden to a minimum) for arbitrary 
measurement models, while also propagating only the first- and second-order statistics of the sensor location. 
Both these requirements are satisfied by an LMMSE estimator. Working with AoA measurements at each receiver, 
the SLE updates are conveniently formulated in the polar coordinates [R 9] centered at that receiver, where R is the 
distance between the source and receiver, and 9 is measured from the x-axis (see Figured] Left Panel). Additionally, 
this choice of polar coordinates makes the LMMSE update optimal under Gaussian measurement models such as 
the LOS model in Section UH 

Each receiver receives the following prior information: SLE fi = [R 9] T and ECM S, where 

^RR ^Re 
^R6 £<9(9 

The AoA estimate at the receiver is 9 with spatial spread gq. Note, henceforth, we use the convention that new 
measurements are represented by 7, prior information by 7, update variables by ?, vectors by . and matrices by bold 
typeface. We desire a linear SLE update, 

fi = p : + K(9-Ap L ), (6) 

where K_ is the Kalman gain and A = [0 1] (only new AoA estimates are available at the receiver). For the 
innovation (9 — Aft) to be orthogonal to the estimate ft, the Kalman gain must be: 

K_ = AY, (a e + At.A T y x . (7) 

Therefore, inserting © into d6), the updated SLE, fi, is obtained as 



R = R + 



See + ere Q 



(8) 



The ECM update can be obtained using © as 



The entries of the updated ECM are given by 

Err = ^RR 



^— = ;=^ + — and Z Re = - (10) 

When the observations 9 have Gaussian errors as in the LOS model (Q3, the LMMSE updates in ([8]) are also the 
optimal minimum mean squared error updates. Under non-Gaus sian error models, the LMMSE updates are the 
optimal linear updates from the perspective of mean squared error. 

The updated SLE jx and ECM S are thus produced in a polar coordinate system centered at the receiver with the 
new measurement. This SLE and ECM must then be transformed to the global cartesian system (a common frame 
of reference) to provide the information in a form accessible to the next receiver. The next receiver then transforms 
these estimates into its own polar coordinate system. We now describe transformations between the local polar 
coordinate and global cartesian coordinate systems. 
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Suppose that the receiver located at X_ c = [x c y c ] T in the global cartesian system computes the new SLE 
pL = [R 9] T with ECM S po ] in polar coordinates. The SLE in the cartesian system, X s = [x s y s ] T , is given by 

2L S = X c + [Rcos(0) Rsm(0)]. (11) 
Therefore, errors in [R 9] T can be mapped to errors in [x s y s } T as 

, * > 

dx a \ _ f cos(#) -£sin(0) \ ( dR 

dy s ) ~ V sm(9) Rcos(9) J \ d9 
and the ECM in the cartesian system is 

Scar = T po i(R,9) Sp i T po i(R,9) T . (12) 

Similarly, for the same receiver at receiving prior SLE X_ s = [x s y s ] T in cartesian coordinates, the prior SLE 
and ECM can be transformed to the receiver's polar coordinates as 

R=\\Kc- 2U I and 9 = arg{X c - XJ, (13) 

and 

Spoi = T CM (R,9) S car T car {R,9) T , (14) 

where T car (i£, 9) = T^(R, 9). Note that the coordinate transformations T po i and T car depend only the measurements 
and SLE at the current receiver and not on the AoA measurement variance. 



B. Sequential Localization Algorithm 

We now describe the steps involved in sequentially aggregating receiver AoA estimates to produce a SLE. Let 
iV receivers be located at Xj. = [x^ yk] T indexed by k, and let the source be at X_ s = [x s y s ] T . Receiver fc's AoA 
estimate is 6%, measured from the x-axis of the global cartesian system. Further, the polar coordinates with receiver 
k at its origin is designated Vk- For ease of exposition, we assume that the receivers are indexed in the order in 
which their estimates are combined. 

The Bootstrap procedure: The AoA estimates of the first two receivers in the combining order, 9\ and #2, are used 
to obtain an initial SLE and ECM to "bootstrap" the Bayesian algorithm (see Figure [T] Left Panel). The estimated 
range of the source from receiver 1 is 

a (2/2 - Vi) cos(0 2 ) - (»2 - x x ) sin(0 2 ) 

K\ = . (ID) 

sin(0i - 2 ) 

Therefore, we can get an initial SLE = [x s y s ] T in cartesian coordinates from fi^ = [R\ 9\} T using (fTTT) . 
Errors in source range estimates from receivers 1 and 2 can be computed from the errors in the AoA estimates as 

dRi \ _ 1 / -Aicob(0i - § 2 ) R2 \ ( d§i 

dR 2 J ~ sin(0i - § 2 ) V ~Ri R-2 cos(^i - 9 2 ) ) V d§ 2 

Using ( fT6l ) and the conditional independence of the estimates 9\ and 9 2 given the source location, the entries of 
the initial ECM Tirr, T<qq, and Srq can be computed in polar coordinates V\. This initial ECM S^] in V\ is 
transformed to the global cartesian coordinates as Scar using (TT2l and fi^. Note that the initial ECM could have 

been computed in V 2 instead of V\, but both would ultimately lead to the same S^. We term these preceding 
steps the 'bootstrap procedure'. 

However, care must be taken in choosing the receivers for the initialization above so as to avoid 'bad' initial 
conditions. It is evident from (031 ) and (fT6l ) that the initial SLE and ECM become unbounded if (#1 — #2) ~ nir, n = 



(16) 
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0, 1, . . ., due to the sin(#i — 62) term in the denominator of both equations. Thus, the error is largest when the 
source is roughly collinear with the two receivers (on the line joining the two receivers), when the nominal AoA 
measurements are close to being parallel or antiparallel. Therefore, the algorithm is always initialized with a pair of 
receivers with 9\ — 62 significantly different from or it, which also ensures that the initial SLE is mostly within 
the ring of receivers. This idea is similar to the concept of dilution of precision in the global positioning system 
||32l that describes the effect of the satellite configuration on the location accuracy. Thereafter, the receivers can be 
combined in any random order and effect of this order on performance is simulated in Section IV-AI 
The sequential algorithm for source localization has the following steps: 

Step 1 (Bootstrap): Estimate initial SLE fi^ and ECM S^j in Vi, using (fT5T ) and (fT6l ). Transform the SLE and 

ECM into the global cartesian coordinates as X^ and £car> respectively, using (TTTT > and ([12). Pass \x} , Scar] as 
a prior to the next receiver. 

Step 2 (Transformation): Let the index of the current receiver be k. Transform the prior 
l \ Sppj ^] in the local polar coordinate system, Vk, using (TT31 and (TBI) . 

Step 3 (Aggregation): Update the prior estimates ] with the AoA estimate of the fcth receiver, §k, 

using the LMMSE procedure in ([8]) and (TTOb . Transform updated estimates [A^,S^] into the global cartesian 
coordinates as 

Step 4 (Termination): If there are unprocessed AoA measurements, pass priors on to the next unaggregated 
receiver and go to Step 2. Otherwise output the SLE and stop. 

Since only the current SLE and ECM need to be passed on from one receiver to the next, this algorithm can 
be implemented in a distributed manner, with each receiver needing to know only its own location and orientation. 
While we consider AoA estimates here, this sequential algorithm is quite general, and can, for example, incorporate 
probabilistic information on the source range obtained from signal strength measurements. Further, the algorithm 
is scalable, in that its complexity grows only linearly in the number of receivers. This scalability is required to 
realize the improvement in localization performance with the number of receivers, details of which are given in 
Section ITIFCl 



C. ML Estimator and Cramer-Rao Bound 

The CRLB is a lower bound on the localization performance of the best minimum variance unbiased estimator. 
For analytical simplicity, we work with the Gaussian AoA model in (Q]), although corresponding bounds for other 
models such as the Laplacian in (0 can be easily derived. For the Gaussian AoA error model, the log-likelihood 
function (within scale factors and constants) for the observed AoA estimates {9i}f =l given the source location 

X = (x,y) is 

L(9 U 9 2 , • • • , N /X) =-J2 (4 "^ 2 ( - ))2 , (17) 

k=l a k 

where 9k (X_) is the true bearing of the source from receiver k and a 1 is the AoA spread (estimation error variance) 
at receiver k. The ML estimator searches for the location X_ that maximizes the log-likelihood function in (TTTT ). 
which is a nonlinear least squares problem. For small o\, the cost function in (fTTT ) can be shown to be approximately 
concave. We have observed numerically that the cost function has an unique global maxima for the parameter values 
of interest and a standard nonlinear least squares solver (such as the "lsqnonlin" function in Matlab® with default 
parameters that is based on [33]) produces the ML estimate. The ML estimate is shown to achieve the CRLB in 
Section IV-Al 

We first construct the Fisher information matrix (FIM) JQQ tnat * s me inverse of the CRLB Ecr(20 = J _1 Q0- 
For the Gaussian LOS model in (fl}, 



^JV sin 2 (6>,Q _V^W cos(6> fc )sin(6> fc ) 
2^k=l alKl 2^k=l 

sr^N cos(6>fc)sin(0 fc ) ^7 

2^k=i z^k=i aim 



T(v\— I ^k=i a jm~ ^k=l ^np k \ n m 

J (^J- 1 v^TV cos(6> fc )sin(6l fc ) cos 2 (6> fc ) > y l °> 
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Y coordinate 00 x coordinate Y coordinate 00 x coordinate 



Fig. 2. Surface plots of the CRLB at various locations on a 500x500 square sensor field surrounded by 8 receivers in a circle for 2° (Left) 
and 4° (Right) AoA spread are shown. The collectors are at (250 + 457cos{27r(fc - l)/8}, 250 + 457 sin{27r(fc - l)/8}), k = 1, 8. 
The CRLB is highest at the center as predicted but some regions close to the receivers also have high bounds. The CRLB contours display 
multiple symmetries due to the equispaced and circular geometry of receivers. 



where 9 = (f?i, 02, . . . , On), and Rf, and 9 k are the range and bearing of the source at location X_ measured from 
receiver k. The total localization error variance is the sum of the variances along the x and y coordinates. The 
lower bound on the total localization error is 

yN 1 

°lm > Tr(£ CR (X)) = — y < 19 > 

Observe that the CRLB is dependent on the location of the source X through {R k }k=v {Ok} N =i (see Figure 0. To 
gain insight into the factors determining the localization error, consider an example of N receivers equally spaced 
on the perimeter of a disc of radius R (as in Figure []]). We simplify the analysis further by computing the bound 
in (fT9l ) for a source at the center of the disc, with the assumption that the spatial spread in the AoA is the same 
for all receivers (we describe a way of handling different AoA spreads at each receiver in Section IIV-Bl i. It can be 
verified numerically that the CRLB is maximized at the center of disc (see Figure and therefore corresponds 
to the worst source position from the standpoint of localization performance. By realizing that 9k = 2irk/N for a 
source at the center of the disc, the lower bound on the localization error at the center of the disc can be computed 

J > NR2a2 ~ 2R2 °* (90) 

" " E^E^sin 2 ^-^^- 2 ^-!)' (20) 

using the fact that sin(0) < 1. Even with a conservative bound, we observe that the localization error reduces at least 
inversely with the number of receivers, thus providing a method of reducing the localization error by increasing the 
number of receivers. The localization error grows linearly with the AoA spread. As seen in Figure [2] although the 
localization error is maximum at the center of the sensor field, the variation of localization error is more complex 
close to the edges of the field. Moreover, AoA measurements are less accurate for far away sources based on (|20j 
and discussion on bootstrapping. Thus, while we consider a ring of receivers for simplicity, the optimal placement 
in order to minimize the worst-case localization error (assuming a fixed area to be covered by a fixed number of 
receivers) is an interesting open issue: we wish to reduce the worst-case distance of a source from the receivers 
nearest it, but we also wish to have a sufficient number of receivers that are near enough to a source. 
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D. Properties of the Sequential Algorithm 

The ML source location is the solution to the non-linear least squares (LS) problem in ([PTl i. The solution to this 
LS problem is asymptotically (in the number of receivers N) consistent and efficient, when the observation noise 
is Gaussian (e.g., LOS AoA error model in (fl3), and the LS solution is asymptotically consistent even when the 
noise is non-Gaussian, but with no guarantees on efficiency [34]. The ML estimate can be computed maintaining 
the asymptotic consistency and efficiency with the following recursive algorithm 11341 : 

v v i i-X/Y ■ s 9L(9 n /X_ n _ l ) o at nn 

An=An-i + J (JU-iJ , n = 2,...,N, (21) 

where X_ n is the source location after combining the first n receiver AoA estimates, J(X n _ 1 ) is the FIM for 
localization with only the first n — 1 receivers, and L(6 n /X_ n _i) is the contribution of the nth observation to the 
log-likelihood function in ( fTTT ). As in Section IIII-CL let [R n 9 n ] T represent the current SLE X n -\ in the polar 
coordinates centered at receiver n. We can now rewrite the recursion for this localization problem as 

1„ =!„-!+ J" 1 (l n -l)( _^coS) )(0n-9n(Xn-l)), U = 2,...,N. 

The above equation has exactly the same form as the LMMSE update in ©, (0, © with the only difference being 
that this is an update in the global cartesian coordinates rather than the local polar coordinates of receiver n. The 
Hessian-like FIM J is also updated at each step and new FIM J(X n ) can be computed using (fT8l) : 

( sin 2 (6>„) -cos(e tl )sin(e„) \ 

-cos(^fs"in(6>„) co^fel) I 

= J(X n -i) + (T-l(Rn,0 n )) H ( J ^T'KRn^n), (22) 



where T po \(R n ,9 n ) is the unitary matrix, in (1121 . used to convert ECM from local polar coordinates to global 
cartesian coordinates. Recalling that the FIM is the inverse of the covariance matrix, observe that the covariance 
update in (l22l is equivalent to the update in the sequential algorithm in © with appropriate transformations to 
cartesian coordinates. Thus, the sequential algorithm is an alternate formulation of this recursive ML computation, 
and therefore inherits the properties of the recursive ML algorithm. However, LS solutions are sensitive to outliers 
in the observations and we present a robust extension of this algorithm next in Section |IV] 



IV. Localization with Non-line-of-sight Channels 

We employ the LOS and NLOS models described in Section [TT] with AOA estimates corresponding to either 
strong multipath components or those that are modeled as outliers, which are chosen uniformly from the set 
of feasible angles. We present an algorithm capable of outlier suppression to handle both the narrowband and 
wideband scenarios. However, we focus on the former (where each receiver produces a single AoA estimate) due 
to its simplicity and later indicate how the algorithm can incorporate multiple AoA estimates at each receiver. 



A. ML Localization under NLOS propagation 

Let a be the fraction of receivers that experience stronger contributions to the received signal from the NLOS 
components than the LOS components. Under the model in dU, the log-likelihood function for the source location 
X , where receiver k has AoA spread a^, is 

N 

£(l/20 = ]T log 

k=l 



a 



a k V27r(l - 2Q(^)) 



exp 



(9 k - e k (x)f 



a 



TT 



(23) 
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where 6 is the vector of observed AoAs. Equation (|23l i can be approximated using log(e a + e b ) 

N 



min(a, b) as 



L(B/X) « - min (4 - e fc (X)) 2 , 6^,* 



(24) 




) 



(25) 



Then the ML SLE X is 



N 



X = argmax L(9/X) ~ argmin 



^min (§ k -e k (X)) 2 ,Ql^ k 



(26) 



k=l 



The ML algorithm in (f26b is a standard LS minimization with angular errors bounded by a threshold max ,fc for 
receiver k. This cost function ensures that there is no incentive to reducing angular errors larger than the threshold 
and therefore, outliers or NLOS estimates only have a limited effect on the SLE. As desired, the ML estimator 
in d26l ) reduces to the ML method for the LOS scenario (see (fTTT l) in the absence of outliers (a — > 0), and the 
threshold G max ^ — > oo. Although the this ML algorithm is computationally prohibitive, we present it to provide a 
heuristic for the outlier suppression algorithm. 

For simplicity of exposition, we describe the outlier suppression algorithm with the same AoA spread for all 
receivers and hence, a common threshold max . To extend the sequential algorithm to the NLOS scenario, we impose 
this constraint, © max , on the largest observed angular error at each step and describe this modified algorithm in the 
following section. 

B. Sequential Aggregation with Outlier Suppression 

In NLOS scenarios, the localization algorithm must identify the largest subset of receivers with LOS channels 
that are mutually consistent, and use these AoA measurements to estimate the source location. However, the 
receivers with NLOS channels (outliers) are unknown and arbitrary in number. This search for the largest mutually 
consistent set of receivers is of exponential complexity, since it must explore every subset of the N receivers. We 
therefore resort to a randomization of the algorithm in Section IIII-B L in which we randomize the choice of the 
first two receivers used in the bootstrap phase. Thereafter, each subsequent receiver's AoA estimate is combined 
ensuring that angular errors in all the aggregated receivers remain below the threshold max from (|25T ) with a k = a. 
By bootstrapping with different pairs of receivers, this randomized algorithm produces source position estimates 
corresponding to different subsets of receivers that mutually agree, leading to a list of possible explanations for the 
observed 9. 

The outlier suppression algorithm is an extension of the sequential algorithm in Section IIII-B I with two key 
differences: First, at each step, the next receiver chosen for aggregation is the one with the smallest angular error. 
The angular error at a receiver is the discrepancy between the bearing of a hypothetical source at the current 
estimated location, X_, and the AoA estimate measured at that receiver, i.e., e(9,X_) = \9 — 9(X_)\. The error e 
can be understood as the empirical estimate of the error in (|26l ). The receiver with the smallest angular error is 
the receiver whose AoA estimate is most consistent with the current SLE. Second, after combining a new AoA 
estimate, the updated source location is retained only if the angular errors for all the aggregated receivers are below 
the threshold max . This ensures that outliers are eliminated using the criterion in (1261 ) and simultaneously, the 
receivers with LOS are used to compute the source location using an LS computation. 

Sequential aggregation with outlier suppression involves repetitions of the following basic steps with multiple 
random bootstraps (for brevity, we use the phrase "combine with 9" to denote the computation of the new SLE 
and ECM using the appropriate coordinate transformations described in Section UlI-AI) : 

Step (Initialization): Set the list of receivers already aggregated A = $ and list of receivers yet to be combined, 



C = {1,... 



N}. 
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Step 1 (Bootstrap): Select a pair of receivers at random, say {i,j} that has not been used previously. Compute 
the initial estimate X_ and ECM S using the bootstrap procedure. Add {i,j} to the list of aggregated (inlying) 
receivers, A = AU {i, j}, and remove it from the list of remaining receivers C = C — {i,j}- 

Step 2 (Angular Error Computation): Compute angular errors e(9,X) over the remaining receivers C, i.e., 

e(e i ,x) = \e i -e i (x)\ vise. 

Step 3 (Candidate Selection and Aggregation): Find the receiver k with the smallest error e(6k,X_)- Combine X_ 
with 9k to obtain the new candidate estimate X! and covariance X'. 

Step 4 (Threshold Verification): If the angular error using the candidate location Jf' is below the threshold for 
all the aggregated receivers, i.e., 

e(^l')<e max , Vl£Au{k}, 

then retain the candidate location and covariance, X_ = X_', and S = S'. Also add receiver k to the list of 
aggregated receivers, A = A U {k}. 

Step 5 (Termination): Remove receiver k from further consideration C = C — {k}. If there are no remaining 
receivers (C = 0) then Stop else goto Step 2. 

Here we assumed that the AoA spread and hence, the threshold is the same for all receivers. As seen from (126V 
the error at each receiver k is only compared against max .fc, the threshold dependent on AoA spread for receiver k. 
Thus, the above NLOS suppression algorithm can be applied to different AoA spreads by replacing the threshold 
in Step 4 by the appropriate maXj fc from (|25l ). 

The above algorithm is repeated with M different random initial conditions in order to detect the source with a 
high probability, and each run produces a likely source location and a confidence (the ECM) in that estimate. In 
a wideband system, if a receiver resolves two arriving paths , the above algorithm can still be used by introducing 
a second virtual receiver at the same location as the original receiver and assigning to it the second arriving path. 
However, it must be ensured that both the original and virtual receivers are not part of the same SLE, since we cannot 
have two different AoA estimates at a given receiver corresponding to an LOS path. This situation is illustrated 
with an example in Section IV-Bi 

Choice of M and max : The performance of the outlier suppression algorithm is determined by the choice of 
the maximum angular error max and the total number of iterations, M. Assuming that the spread in AoA is known 
at each receiver, the threshold max can be computed using d25l) . if the fraction of receivers with strong NLOS 
components a is known. In practice, the largest expected fraction of receivers with strong NLOS components, a max , 
can be set based on knowledge of the propagation environment and the worst case NLOS scenario under which we 
wish to operate. This value of max obtained using (l25l) is conservative for lower levels of multipath scattering, as 
max is monotonically decreasing with a. While choosing a smaller max might prevent some "good" LOS AoA 
estimates from being utilized, it also ensures that NLOS measurements do not corrupt the SLE. 

During the bootstrap phase, two receivers are chosen randomly to seed the sequential algorithm. Success in the 
sequential estimation depends on selecting two receivers with LOS channels to initiate the algorithm. We have 
observed from simulations that the algorithm always converges to a solution in the 'vicinity' of the bootstrap 
location, and therefore an estimate using almost all the LOS estimates is produced, if the algorithm is bootstrapped 
with two receivers having LOS channels. Hence, we hypothesize that bootstrap failure is the predominant cause of 
localization failure (we verify this numerically in Section IV-BI ). We now try to estimate this probability of failure 
as a function of M, the number of iterations of the outlier suppression algorithm. Suppose the outlier suppression 
algorithm is repeated with M different seeds, the probability of failure in the bootstrap phase is the probability 
that at least one receiver is an outlier in each of the M attempts. The total number of bootstrap pairs is P = ( 2 ) 
and number of pairs with at least one NLOS receiver is K = P- ( L(1 "^ )7VJ ), where [(1 - a)N\ is the number of 
receivers with LOS channels. Then, 

{LA if M < K 
(L) ~ (27) 

if M > K 
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When the receivers resolve multiple arriving paths, the above probability of failure computation is modified by 
replacing the number of receivers by the total number of resolved paths. Depending on the probability of failure 
acceptable in the system, d2Tb is used to choose the number of randomizations of the algorithm that are necessary. 

A practical issue of interest is the choice of the number of randomizations M for different total number of 
receivers N to achieve a fixed probability of bootstrap failure under similar NLOS propagation environments (fixed 
a). Rearranging the tight upper bound on the probability of bootstrap failure, presented in the Appendix as (130) , 
we get an tight upper bound on the possible value of M as 

^ log(P(bootstrap failure)) 

M — i — n n \9\ ' 

log(l - (1 - a) 2 ) 



using the fact that for large N 



— = i - A ? )_ x _ (1 _ y. 



We observe that the probability of bootstrap failure is independent of N. Thus, the outlier suppression algorithm 
with complexity 0(MN 2 ) is still only quadratic in the number of measurements. 



V. Numerical Results 

We study the performance of the proposed algorithms via Monte-Carlo simulations. The simulation setup is 
as follows: We consider a circular field of unit radius with N equally spaced receivers along the perimeter. As 
seen from our analysis in (|20l , the localization error grows linearly with distance from the receiver. Therefore, by 
selecting a field of unit radius, we obtain scale-invariant (only dependent on dimensionless quantities) measures of 
performance. Let the receivers be located at [cos(2-7r(A; — 1)/N) s'm(2ir(k — 1)/N)] T , k = 1, . . . , N. Each receiver 
measures the AoA of the signal received at its antenna array and representative AoA estimates are generated 
according to the models described in Section |II] We assume throughout that all the receivers have the same AoA 
spread for convenience, although the algorithm does not require this, and that only a single source transmits at any 
given time. 

The AoA spreads under different environments reported in literature are listed in Table U The effective AoA 
spread for our NLOS models in © is the quantity reported in most studies. This effective spread is greater than the 
AoA spread a for the Gaussian LOS model alone due to the mixture with the uniform distribution corresponding 
to LOS blockage. Therefore, accounting for the available information on propagation environments in Table U we 
consider AoA spreads a in the range 1 — 10° for numerical studies of our algorithm. We use the CRLB (TT9T l as 
a performance benchmark, but this bound is dependent on the true location of the source. In order to have a fair 
comparison, we run equal number of iterations on each of 25 candidate source locations, and compare the total 
rms localization error against the average of the CRLB at those 25 locations. 



TABLE I 

AOA SPREADS FROM REPORTED MEASUREMENTS IN DIFFERENT ENVIRONMENTS 



Environment 


AoA spread 


Reference 


Outdoor - Urban 


~ 10° (LOS + NLOS) 


Figures 2 and 3, Klein et al. [35] 


Outdoor - Urban 


~ 15° (LOS + NLOS) 


Figure 5, Thomas et al. |36| 


Indoor - LOS blockage 


20 - 30° (NLOS) 


Spencer et al. [37| 


Indoor - WIMAX 


10° (LOS), > 50° (NLOS) 


Figures 6 and 7, Akdhar et al. [38] 


Outdoor - Urban 


3 — 20° (LOS) depending on distance 


Figure 10, Chen and Asplund [39| 



A. Performance under LOS scenarios 

The performance of the sequential algorithm in Section IIII-BI for N = 8 receivers is shown in Figure [3] The 
algorithm achieves the CRLB for small angular estimation errors (in simulations for as few as 6 receivers), but 
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the performance deteriorates for large spatial spreads in AoA. The optimal ML estimator, in Section UlI-CI can be 
shown to be approximately convex, and the likelihood function has a unique maxima. The sequential algorithm, 
which is an instance of a stochastic approximation algorithm [34], has a tendency to get stuck in local minima that 
arise with large AoA spreads; at large spreads, there are multiple locations where subsets of estimated AoAs are 
(roughly) concurrent, but no single location close to all the directions is evident. 

This problem is solved by selecting the most likely SLE (i.e., has the highest ML cost in ( fTTT )) from multiple runs 
of the sequential estimation using different pairs of receivers to bootstrap. From the simulations in this section, it 
appears that the choice of a random combining order for the AoA measurements does not incur any significant loss 
in performance against the CRLB. In reality, the nonlinear coordinate transformations at each step in the sequential 
algorithm are dependent on the current SLE making the final estimate weakly dependent on the specific order. 
Nevertheless, the remaining gap to the ML performance can be closed by using multiple random bootstraps and 
selecting the estimate with smallest ML cost, as shown in Figure [3] and conforms well with the randomization 
framework in the NLOS suppression algorithm in Section HV-Bi 




i , , , , , , , 1 

2 4 6 8 10 12 14 16 18 
Std. deviation of AoA errors (in degrees) 

Fig. 3. Performance of sequential algorithm under LOS setting for N = 8 receivers is compared against the CRLB (solid) and the ML 
estimate (dashed). The performance improvement using multiple random bootstraps (1,2 or 3) is also shown. 

Incorporating range information: The localization resolution can also be improved by using additional modalities 
such as received signal strength (RSS)-based ranging. In order to demonstrate the flexibility of our algorithm to 
also leverage range information, we compare our system, utilizing only AoA measurements, to one where all the 
receivers additionally obtained RSS-based range information. The RSS-based range estimate, R, is modeled as a 
log-normal random variable (independent of the AoA estimate) as suggested in (7) within fiOl : R = i?exp ( ^g^ ), 
where R is the true source range, r is zero-mean unit variance Gaussian random variable, and /3 is the path 
loss exponent. For the simulations, we chose /3 = 3 and the variance of the random variable R was taken to be 
the measurement error variance. The same update formulation in Section IIII-AI can be extended to include range 
information as well. Given the new range and AoA estimates jl = [R 9] T with covariance S, the Kalman gain K 
is computed as K = S(S + and thus, the new SLE is /t = Jl + K(/i — Jl). The updated covariance matrix 

is 5T 1 = IT 1 + XT 1 . 

In Figure 01 we show the normalized localization performance for three scenarios where the percentage errors 
in the AoA measurement are below (bottom curves), approximately equal (middle curves) and greater (top curves) 
than the percentage errors in the range estimates, respectively. As expected, as the quality of the range information 
improves, the gains from using range information are greater. Further, when the range measurements are less 
accurate than the AoA estimates, performance with hybrid range and AoA information can be approached using 
AoA measurements only, with moderate increases in the the number of receivers, N. This is possible as location 
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AoA + Range (a = 2 ) 

AoA only (a = 2°) 

'+- AoA + Range (o = 4°) 
■> — AoA only (a = 4°) 
O- AoA + Range (cs = 8°) 
«— AoA only (o = 8°) 




10 12 
No. of receivers N 



Fig. 4. Comparison of localization performance of the sequential algorithm in a system with only AoA information and one where all 
nodes have RSS-based range estimates in addition to AoA measurements for different numbers of receivers. The receiver range estimate is 
modeled as a log-normal random variable parameterized by SNR = 12 dB. 



errors due to both RSS-based range estimates and AoA estimates scale linearly with the range of the source from 
the receiver, allowing the substitution of range information by AoA estimates from additional receivers. Moreover, 
this similarity in scaling makes the improvement in localization resolution with number of receivers independent of 
the modality used (see Figure @). Finally, a key attribute that we reiterate here is the capability of this framework to 
easily incorporate other sensing modalities. Although we do not explicitly demonstrate here, outliers in the range 
estimates can be handled similarly in the algorithm in Section IIV-BI by adding a range threshold. 

B. Performance under NLOS scenarios 

In this section, we explore the capabilities of the outlier suppression algorithm under the narrowband and wideband 
multipath models. 




3 4 5 6 7 
AoA error variance (in degrees) 

Fig. 5. Localization performance with NLOS suppression with TV = 8 receivers for different fractions, a, of receivers with NLOS 
channels with optimal thresholds chosen according to d25t . The location error variance (solid line) is compared against the ML error variance 
(dash-dotted line). 
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1 ) Effect of multipath in a narrowband system: We first present numerical results for the narrowband multip ath 
model, where each receiver can only resolve paths spatially and each receiver produces a single AoA estimate 
corresponding to either the LOS path, or the reflected and scattered multipath. In Figure |5J we compare the 
localization error using the outlier suppression algorithm and the optimal ML estimator in Section ITV-AI for a system 
with 8 receivers. For different fractions of outliers, a, the outlier suppression algorithm is run with the threshold 
©max chosen according to ( |25T ), while the number of random bootstraps is chosen to ensure that the probability of 
bootstrap failure is less than 1CP 3 . This resulted in a choice of the number of random seeds, M = {4, 7, 11, 15}, 
for a = {0.125,0.25,0.375,0.5} (or {1,2,3,4} outliers) using (|28l ). The algorithm puts out multiple solutions, 
one corresponding to each random initialization. After pruning out the estimates that placed the source outside the 
circular field, the ML cost function in (|23l l is used to select the most likely estimate. The ML estimate was obtained 
by brute force minimization of the same cost function. 

The algorithm performs very close to the optimal ML estimator for the entire range of AoA spreads. However, 
it is interesting to note in Figure [6] that the ML estimate does perform significantly worse compared to the ML 
estimate using only the good LOS AoA estimates. This additional loss is the cost of identifying the NLOS receivers 
and as the angular spread increases, it becomes progressively more difficult to differentiate between the LOS and 
NLOS AoA estimates. 
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1 outlier (K) 
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2 outliers (K) 

3 outliers (U) 
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Std. deviation of AoA errors (in degrees) 

Fig. 6. The ML localization performance for an N = 8 receiver system with different numbers of receivers with NLOS channels for the 
two situations when the receivers with NLOS channels are known (K) or unknown (U). The gap between the performance when the NLOS 
receivers are known (dash-dotted line) and are unknown (solid line) represents the performance penalty for outlier identification. 




In practice, since the fraction of outlying receivers, a, is unknown, the algorithm is operated with a threshold 
chosen using an upper bound on this fraction, which in our simulations is a max = 0.5. In Figure |7J for different 
numbers ({0, 1, 2, 3, 4}) of NLOS receivers, the outlier suppression algorithm achieves close to ML performance for 
smaller actual fractions of outliers even for this conservative choice of threshold. Thus, we can conclude that this 
approach is quite insensitive to the exact choice of threshold, max , which adds to the robustness of this approach. 
The simulation results also indicate that the algorithm in Section HV-B I is approximately ML (AML). 

In Figure [8j the observed probability of failure and the expected probability of failure from (|27l i of the outlier 
suppression algorithm are plotted for N = 8 receivers. A failure is declared if the SLE is outside a circle of 
radius three times the standard deviation of the ML algorithm. When the SLE errors are normally distributed 
N(0, u 2 ), as is expected from our analysis of the sequential algorithm in Section Hll-D[ the rms localization error 
is Rayleigh(^). Thus, the probability of AoA measurement "noise" alone causing the estimate to lie outside the 3v 
circle is exp (— (3v) 2 /v 2 ) rj 10~ 4 . Hence, failures due to bad bootstraps are significant only when the observed 
failure probability is of the order of 10~ 3 and above. We observe from the figure that the expected failure probability 
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0.45 




123456789 
AoA error variance (in degrees) 

Fig. 7. Localization performance with NLOS suppression for an TV = 8 system in a realistic scenario with threshold chosen with the 
maximum expected outlier fraction, a max = 0.5 for different fractions, a, of receivers with NLOS channels. The location error variance 
(solid line) is compared against the ML error variance (dash-dotted line). 



is greater than the observed probability over almost the entire range of interest, but plateaus around 3 x for 
all three fractions of outliers. This, we believe, is due to the fact that localization errors are not strictly normal in 
the presence of outliers leading to a slightly higher failure rate due to AoA measurement "noise". We can safely 
conclude, therefore, that bootstrap failures dominate the failure probability over the entire range of interest. 
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Fig. 8. The comparison between the observed probability of failure and the expected probability due to bootstrap failures alone for a system 
with N — 8 receivers with different fractions of outliers a. The algorithm is defined to have "failed" if the SLE lies farther than three times 
the standard deviation of the ML estimator from the true location. 



2) Effect of alternate narrowband multipath models: Although our algorithm is AML under the LOS and 
narrowband multipath models developed in Section HH it is of interest to examine the sensitivity of the NLOS 
suppression algorithm to these models. To this end, the algorithm was simulated, shown in Figure |9] with two 
heavy-tailed AoA error models, namely the Laplacian (in (|2}) and Cauchy (the standard deviation corresponds to 
the shape parameter here), for the AoA estimation errors with the nominal parameters values in Figure [7] we did 
not compute the ML solution for the Cauchy model as it is mathematically intractable. Under the Laplacian model, 
which has smaller tails, the algorithm attained the optimal ML performance for that model. But, with a Cauchy 
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model, the heavier tail generates many more outlying AoA estimates leading to larger estimation errors, and the 
observed performance is equivalent to that of our nominal model in (@]) for a = 0.375. 




123456789 
Std. deviation of AoA errors (in degrees) 



Fig. 9. Localization performance of NLOS suppression when AoA estimates are obtained from Laplacian and Cauchy models with threshold 
chosen with the maximum expected outlier fraction, a m ax = 0.5, as before. 

3) Ray-tracing Illustrations: We now compare the performance of the outlier suppression algorithm in the 
narrowband and wideband setting with the following two examples. In both examples, we use a virtual point source 
(shown as a solid circle) model to trace multipath generated by the reflectors; the LOS path from the virtual source 
to a receiver corresponds to the nominal direction of arrival of the reflected signal from the true source. In all 
examples, AoA estimation errors of standard deviation 0.5° are added to both the nominal direction of arrival of 
the LOS path and the multipath from the virtual source. 

Narrowband multipath setting: In Figure \\0\ four receivers (squares) attempt to locate source A ('+' sign) in 
the presence of two reflectors (solid gray lines). One wall blocks the LOS path to receiver R4, causing R4 to only 
receive multipath reflected by the second wall. The receivers Rl and R2 have LOS AoA estimates, while receiver 
R3 receives both the LOS path and multipath from the virtual source (i.e., reflecting wall). The narrowband receivers 
generate one AoA estimate each, corresponding to the superposition of all the arriving paths. Thus, receivers Rl 
and R2 have very reliable LOS AoA estimates, R3 estimates an AoA with a large spatial spread and receiver R4 
'sees' an outlying AoA measurement. The SLEs from multiple runs (crosses) of the outlier suppression algorithm 
is shown in Figure [10] The availability of reliable LOS estimates from Rl and R2 produces good estimates of the 
source location by eliminating the outlying estimate from R4 and also prevents the algorithm from mistaking the 
virtual source to be a real source. Localization performance at three other locations, B,C and D, is also shown in 
Figure [TO] Good performance can, therefore, be expected in the narrowband setting even in the presence of LOS 
blockage if there are sufficiently many receivers with LOS paths to the source. However, there are situations, like 
when all the receivers experience NLOS propagation, where the narrowband system performs poorly. We elaborate 
on this issue in the following example. 

Wideband multipath setting: As described in Section [TTJ when the source transmits a wideband signal, the receiver 
receivers can additionally resolve arriving paths in time, since reflected multipath components suffer delay with 
respect to the LOS path. This leads to multiple AoA estimates at the receiver corresponding to different directions 
and times of arrival. However, we cannot always reliably conclude that the earliest arriving path is LOS. Instead, we 
choose to apply our outlier suppression algorithm to all the estimated AoAs and allow the algorithm to eliminate 
NLOS AoA estimates as outliers. We illustrate this with an example in Figure [TT] A source ('+' sign) is situated 
between four receivers (squares) and a wall (solid gray line). Under the wideband setting, in Figure [TT] each receiver 
generates two AoA estimates, one due to the LOS path and the other due to the reflected NLOS path, while with 
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Fig. 10. The multipath suppression capability of the algorithm in a narrowband setting is illustrated with an example. With the source at 
location A, receiver R4 experiences LOS blockage and receives only reflected multipath from the source. The output of multiple runs of the 
outlier suppression algorithm with M = 7 is shown for location A only. Location B uses only 3 receivers as R4 does not receive any signal 
from the source, location C uses all the receivers with LOS and location D has similar ray traces to location A. The rms error at these four 
different locations under this configuration are 18.42 at A, 2.7 at B, 2.57 at C and 14.22 at D (individual estimated locations not shown for 
B, C and D). 

the narrowband setting, the receivers estimate the AoA as a power-weighted superposition of the two directions of 
arrival. 

The output of multiple runs of the outlier suppression algorithm is plotted for the wideband and narrowband 
system. In the wideband system, the virtual source location is identified as a likely position in addition to the true 
source location and the NLOS algorithm is not capable of differentiating between true sources and virtual sources 
arising due to correlated multipath. But in practice, the knowledge of the environment can be used to eliminate 
infeasible estimates such as source locations behind the wall. On the contrary, in the narrowband scenario, the source 
is located in the region between the true and the virtual sources, as the multipath, in effect, increases the spatial 
spread in the LOS AoA estimates and degrades performance greatly. It is apparent that the capability to resolve 
multipath is essential to achieving satisfactory performance in NLOS environments and helps on two counts: First, 
resolving multiple incoming paths reduces the effective spatial spreading on each path. Second, multiple estimates 
at each receiver increase the total number of "good" LOS measurements to estimate the source. 

VI. Conclusions 

The sequential algorithm nearly achieves the CRLB in an LOS environment, and forms the building block for 
our outlier suppression algorithm for NLOS environments. Our algorithms have at most quadratic complexity in 
the number of measurements, are amenable to distributed implementation and can be easily modified to incorporate 
other sensing modalities such as RSS-based range measurements. For LOS environments, the scaling behavior of the 
localization error variance is linear in the AoA spatial spread (which determines the variance of the AoA estimates) 
and the size of the coverage area, and is inversely proportional to the number of receivers. Since we present results 
in terms of scale-invariant quantities, we can estimate the localization accuracy that can be obtained under any 
scenario as (rms error from simulations ) VArea covered/7r. Thus, for instance, in an outdoor environment of area 1 
km 2 with intermediate level of local scattering (we assume 4° spatial spread from Table H), the obtained localization 
accuracy using 8 receivers is about 15 m using simulation results in Figure |4l On the other hand, with an angular 
spread of 9° in the LOS path and 10% outliers in an indoor WIMAX environment [38], an accuracy of 2 m can be 
obtained in a room of size 50x50 m using 8 receivers (see Figure [5]>. In NLOS multipath scenarios, although the 
outlier suppression algorithm is approximately ML, its performance is heavily dependent on the specific type of 
environment and the capability of the receivers to resolve the contributions of the LOS path and NLOS multipath in 
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Fig. 11. The working of the localization algorithm in the presence of a virtual source due to a perfect reflector (wall) with wideband and 
narrowband source transmissions. The output of multiple runs of the algorithm is shown. 

the received signal. In a narrowband system, where each receiver only resolves a superposition of the arriving paths, 
the algorithm can suppress outlying AoA estimates, if there are sufficient number of receivers with reliable LOS 
estimates. However, the capacity of a wideband system to estimate AoA from the LOS and multipath individually is 
vital in settings where all the receivers experience NLOS propagation. There is a localization performance penalty 
for having to "find" the outliers, which becomes progressively worse as the fraction of outliers in the measurements 
increases. 

Broadly speaking, the key issue for further investigation is to relate physical models of the transceivers and the 
propagation environment to the AoA estimation models that form the basis for our algorithms and analysis. This 
includes understanding the effect of the antenna array geometry, the placement of the receivers, and the propagation 
environment. In particular, for canonical environments of interest (e.g. flat outdoor terrain, urban environments, 
indoor office environments, warehouses), it is of interest to develop models, given a deployment of receivers, for 
SNR variations at the receivers (due to variations in their distances from the source and multipath fading) and the 
fraction of receivers with blocked paths. This in turn would provide a framework for optimizing the deployment 
of receivers. 
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Appendix 

In this appendix, we derive the bounds on the probability of the bootstrap failure in (127T ). As defined earlier, 
P = ( 2 ) is the total number of bootstrap pairs and K = P - ( L(1 ~2 )JVJ ) is the number of pairs with at least one 
NLOS receiver, where [(1 — a)N\ is the number of receivers with LOS channels and a is the fraction of receivers 
with NLOS channels. 



P (bootstrap failure) = if M < K 

( K-P\ ( K-P\ ( K-P \ 
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In order to obtain an upper bound on the bootstrap failure probability, we replace each ratio of the form by 
a larger fraction K p P (note that K < P by definition): 

M 

P 



^(bootstrap failure) < ( — ) . (30) 



Similarly, replacing the ratios p_f by a smaller fraction p^/j^ in d29l ), we get a lower bound, 



JU+l 

K - M + l x " 



^(bootstrap failure) > I — — - J . (31) 

The lower and upper bounds clearly converge if M <C K < P, which occurs when N is large. Moreover, these 
bounds hold only for M < K, as the probability of bootstrap failure is zero if M > K. 
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